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It is discussed the decoherence problems in ensemble large-scale sòlid state NMR quantum com- 
puter based on the array of 31 P donor àtoms having nuclear spin / = 1/2. It is considered here, as 
main mechanisms of decoherence for low temperature (< 0.1 K), the adiabatic processes of random 
modulation of qubit resonance frequency determined by secular part of nuclear spin interaction with 
electron spin of the bàsic àtoms, with impurity paramagnetic àtoms and also with nuclear spins of 
impurity diamagnetic àtoms. It was made estimations of allowed concentrations of magnètic im- 
purities and of spin temperature whereby the required decoherence suppression is obtained. It is 
discussed the random phase error suppression in the ensemble quantum register bàsic states. 
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Introduction 



There are five bàsic criteria for realization of a large-scale quantum computer, which can outperform all tra- 
ditional classical computers 0. Let us take a brief look at these criteria in respect to the large-scale NMR 
quantum computers. 

1. For any physical system, which presents large-scale quantum register, it is required of a sufficiently large 
qubit numbcr (L > 10 3 ) and each qubit should be separately identifiable. 

One such example of this register is solid-state homonuclear system, in which nuclear spin containing identical 
àtoms are housed at regular intervals in a natural or an artificial solid-state structure. 

2. There is a need to provide the conditions for preparation of initial and bàsic quantum register state. 
For a many-qubit NMR quantum computer the quantum register state initializing can be obtained by going to 
extra- low nuclear spin temperature (< lmK at fields of order of several Tesla). 

3. The decoherence time of qubit states Td should be at least up to 10 4 time longer than the "clock time". 
It has value of order of several seconds for NMR quantum computers. The decoherence suppression is one of 
the important problems in realization of a large-scale quantum computers. 

4. There is a need to provide during a decoherence time the implementing of a set of quantum lògic operations 
determined by a lògic unitary transformation. This set should contain certain set of the one-qubit and two-qubits 
operations are shielded from random errors. The electromagnètic pulses that control the quantum operation 
should be performed with an accuracy of better than 10 -4 -10 -5 . 

5. There is a need to provide the accurate and sensitive read-out projective measuremcnts of individual 
qubit states. This is another of the important and the most hard problems. 

One of the pioneering schcmes of large-scale NMR quantum computer with individual access to qubits was 
proposed by B.Kane H |J. 

There are four peculiar difhculties in implementing of such quantum computer: 

1. First of all, signal from the spin of an individual atom is very small and it is required of a highly sensitive 
single-spin measurements. 

2. For initialization of nuclear spin states it is required to use very low nuclear spin temperature (~ mK). 

3. It is required to use the regular donors and gates arrangement with high precision in nanometer scale. 

4. It is necessary to suppress the decoherence mechanism defined by thermal fluctuations of gate voltage 
(this mechanism will be denoted as external). 

As an alternative, we have proposed the variant of an ensemble silicon-based quantum computer H, [|, |(|. 
One would expect that with the ensemble approach when many independent "molècules" of Kane's type, like 
in the bulk-ensemblc liquid quantum computers, work simultaneously, the qubit state measurements seem to be 
greatly simplified. Furthermore, it seems to be solved certain problems associated with possible noncontrollcd 
continuous random phases of individual qubits. 



This article discusses only the decoherence mechanisms of one qubit quantum states in ensemble quantum 
register based on the regular array of 31 P àtoms having nuclear spin 1=1/2. It is considered, as main 
mechanisms of decoherence for the very low temperature, the adiabatic processes of random modulation of 
qubit resonance frequency without spin flips produced by fluctuating local magnètic fields which are determined 
by secular parts of interactions of nuclear spins with electron spin of the bàsic àtoms, with impurity paramagnctic 
àtoms and also with nuclear spins of impurity diamagnetic àtoms. We have named this mechanisms as infernal. 
It was made estimations of allowed concentrations of magnètic impurities that are necessary for decoherence 
suppression. It is discussed also the random phase error suppression in the initializing proccss for ensemble 
quantum register states. 

1 Mechanisms of qubit state decoherence in NMR solid-state quan- 
tum registers 

The full electron spin polarization will be practically achieved at 

ls hB/kT » 1, (1) 

where 73 = 176radGHz/T is giromagnetic ratio for electron spin, h = 1.05 • 10 -34 Js/rad, B is induction of 
external magnètic field, T is temperature of environment (of electron spins). Thus for P< = 2T, T = 0.1K 
(B/T = 20 T/K) we obtain j s hB/kT = 27. 

The electron and nuclear longitudinal relaxation processes, that are followed by spin flips in four energy 
level system of phosphorus àtoms doped in silicon, have been previously investigated experimentally in Q]. 
The electron longitudinal relaxation times at low temperatures were found to be extremely long (hours) and 
independent of phosphorus concentration below ~ 10 16 cm~ 3 . The nuclear longitudinal relaxation time Tji were 
found to be equal to 10 hours 0. 

The extremely long longitudinal relaxation times of the electron and nuclear spins imply the possibility to 
produce a long-lived nuclear nonequilibrium initialized state for lattice temperatures (say, for T = 0.1 K) when 
only the electron spins are polarized. The required initializing of nuclear quantum states (near-full nuclear 
polarizations) can be obtained by a short duration deep cooling to Tj < 1 mK only for nuclear spin system 
without deep cooling of the whole sample. There is the possibility to achieve it by indirect cooling of nuclear 
spin system using dynamic nuclear spin polarization techniques (such as the known Abragam's solid-state effect) 

i- 

The relaxation of nonequilibrium state of the nuclear spin system represented by the product of independent 
(nonentangled) one-qubit states, owing to the interaction with isotropic environment, shows two processes. One 
is a slow establishmcnt of equilibrium state associated with dissipation of energy. For it the diagonal elements 
of density matrix decay with characteristic longitudinal (spin-lattice) relaxation time Tu. The decay of non- 
diagonal matrix elements called decoherence of quantum states is characterized by a decoherence time Td or 
transverse (spin-spin) relaxation time T±. The longitudinal relaxation times Ty in the case of nuclear spin of 
31 P àtoms as qubits is defined mainly by thermal modulation of hyperfme interaction accompanied by spin flips. 
It is usual that for sòlids Tj_ <C Tji . 

We will not consider here the external decoherence process due to gate voltage noise. It was made in [0, 0, 0], 
The infernal adiabatic decoherence mechanisms due to a random modulation of qubit resonance frequency, 
produced by local fluctuating magnètic fields without spin flips, seem to be the leading. 

2 Semiclassical model of adiabatic decoherence of one-qubit state 

We will consider a long-lived non-equilibrium qubit state when diagonal elements of density matrix may be 
treated as a constant. 

The random modulation of resonance frequency Ao>(í) that causes the dephasing of a qubit state are deter- 
mined by the random phase shifts 

ip(t) = [ Au(t)dt. (2) 







The one-qubit density matrix of pure state in rotating frame with non perturbed resonance circular frequency 
will be 



p(t) = 1/2 



1 + P Z P_exp(i<p(í)) 
P+ exp(-íp(t)) 1-P, 



(3) 



p, 2 = l. 



where P± = P x ± iP y , P x , P y , P z are Bloch vector components of length P - 

By treating the resonance frequency modulation as Gaussian random process after averaging (|3|) over phase 
distribution with (ip(t)) — we obtain 



<P(*)> = 1/2 



1 + P z P_ exp(-r(t)) 

P+exp(-r(í)) l-P, 



where 



r(t) = 1/2 • ^(jf Au(í)dí) 2 ^ = jf (i - t) (Aw(r)Ac;(0)> dr, 



(4) 



(5) 



/(i) = (Aw(í)Acj(O)) is the frequency correlation function of a random process, which is characterized by 
variance (Au;(0) 2 ) and correlation time tq such that for t > tq (Acj(í)Aw(O)) 0. For T(t) > the averaged 
density matrix presents a mixed quantum state with two non-zero eigen states 



1/2 • 1 ± Jl (P? + P y 2)(l - exp(-2r(í)) 



(6) 

and the populations of states p± = 1/2(1 ± P z (0)) at T(t) => oo. 

Thus, the adiabatic decoherence problem is reduced to the determination of the function T(t) or the corre- 
lation function of random frequency modulation. 

In the case of an ensemble quantum register there is a need to average the one-qubit density matrix and 
correlation function over ensemble of independed equivalent spins-qubits. 



3 The nuclear spin states decoherence due to hyperfine interaction 
of nuclear and electron spins 

In this case the modulation of nuclear spin resonance frequency Au; (i) is determined by the secular part of 
hyperfine interaction: 

Au(t) = A{t)S z (t) - A (S z ) . (7) 

The most studied exactly solvable one-boson models of adiabatic decoherence that is not followed by spin flips 
are not adequate |^0| . The decoherence that is described by the direct one-phonon process with spin flips 
is non-adiabatic slow relaxation process. It is characterized by a rate comparable to 1/Tji- To describe the 
dcphasing or pure adiabatic decoherence we will take in to account one of the Raman two-phonon processes, 
which describes the hyperfine interaction constant modulation A(t) due to the phonon scattering without of 
spin flip. In interaction representation we will have 

A(t) = A [l + y^g ml b+b] exp(i(u; m - wi)t) + g* nl bfb m exp(-i(w m - U)i)t)], (8) 

m#l 

where Aq — 725radMHz is hyperfine interaction constant for àtoms 31 P at T = 0, òj + and ò m are the emission 
and absorption of phonon operators with freqüències u)\ and u m , g m \ ~ the coupling 

constants of hyperfine parameter with acoustic phonons, dimcnsionless parameter £ in the worse case is ~ 1, v 
is the velocity of sound, N is the total number of àtoms in crystal (for simplicity we assume the simple cúbic 
lattice), M is the atom mass. 

From (Q) we obtain the sum of two independent terms, which correspond to the two mechanisms of adiabatic 
decoherence due to hyperfine interaction: 

Aw(í) = Aw s (í) + Auj h (t) = A (S z (t) - (S z )) + 

+1/2 • A o y^g m iòi6+ exp(iu; ml í) + g^b+b m exp(-iu> m \t), u m i = uj m - u\. (9) 

At first, let us consider the first term in (|^). The correlation function is determined by the fluctuations of electron 
spin polarization and depends on electron resonance frequency wg, longitudinal n (hours) and transverse Ti 
relaxation times. In adiabatic case ws = > 1/ri 3> l/Ta and we will obtain: 



(A s o;(í)Aws(0)) = (Au;i) • exp(-í/n), 



(10) 



where 

(Au| ) = AK(Sl) - (S z ) 2 = A 2 (l - tanh 2 ( 7s nB/fcT))/4. (11) 

Now 

r(t) = (Aüj 2 s ) 7? (t/n - 1 + exp(-í/ri)). (12) 

For n » 10 4 s and í - T d = 1 s, 1 < (Aw§) ?f < (n/T d ) 2 we have the non-Markovian random process (slow 
damping fluctuations). In this case 

r(í) = (A^)í 2 /2, (Ac 2 )«i/t 2 . (13) 

Let us write the requirement for the decoherence time for "/shB/kT 3> 1 in the form 

1/T 2 « A 2 Q (1 - tanh 2 ( 7s nS/fcr))/4 « A 2 exp(- 7s ?LB//cT) < ls~ 2 , (14) 

from which we find that the required decoherence suppression will be achieved at sufficiently large ratio B/T > 
30 T/K. For B/T = 20 T/K we have T d - 10~ 3 s. 

For decoherence description due to the second term in (|^) it is convenient to use directly the function T(t). 
With 

f A^ b (í)dí - 1/2 ■ Ao^iU^b+h + CmW^m), 
J ° m,l 

Cmi(í) = g m ii(l - exp(zu; m ií)/cj m i, (b£fim) = n(w m , T) = (exp(cj m /T - 1) _1 

we will obtain 

r(í) = 1/4- ^jf Aw b (t)dt) ^ - 

= 1/4 • A 2 V(n/2A/« 2 7V) 2 • (w mWl )(2n(w m) T)n(cü h T) + n(cü m , T) + n(oj h T)) • g^J^) , (15) 
Let us go now from sums to integrals: 

g]\j r^ 1 

X/"-^ 03 / ■•• Ci; m dw m, (16) 
m ^° 

where Í2 = v(6t:) 1 ^ 2 /a — k<d/h is Debye frequency, a is the lattice constant. 
Taking into account that for uj m \ => (near-elastic scattering) 

sin 2 (o; ml t) 

2 ~ 7r|í|ó(u; m i), (17) 

we now obtain 

r(í) = ^A 2 (^/M V 2 ) 2 (T/e) 7 (fce/n) £ /T ^^ dx \t\ = \t\/ Td . (is) 

For silicon: 6 = 625 K, a = 5.4 • 10~ 8 cm, v « 5 • 10 5 cm/s, Aí = 0.46 • 10- 29 Js/cm 2 . 
At low temperature T/Q < 10~ 3 for T d estimation we obtain: 

i/T d » ^A^/Mi; 2 ) 2 (fce/?i)(r/e) 7 - 3/4 • io 4 (T/e) 7 s - 1 < 1 s _1 . (19) 

8 

Thus, we see, that at low temperatures the phonon mechanism have an insignificant effect as compared with 
the mechanism of electron spin fluctuations. 



4 The nuclear spin states decoherence due to interaction with mag- 
nètic impurity àtoms. 

An alternative reason for the modulation of individual nuclear spin resonance frequency is the secular part of 
their dipole-dipole interaction with random distributed paramagnetic center (impurity atom and defect) in the 
substrate, which for the magnètic impurity àtoms in S-state is of the form: 

Acj(t) = 7I B Z (Í) = £ • 7l7s ?i • ^ 1 ~ 3 r 3 i2/r ' 2 ■ (S z (r u t) - <S z (rO)), (20) 

i 1 

where /io = 0.47rT 2 cm 3 /J, T{ is radius-vector of the impurity atom distance. 
In this case, once again, we have 

(Aw(i)Aw(Q)) = (Aw 2 ) • e X p(-í/T 1)inip ), (21) 

where for non-correlated homogeneous spatial distribution of impurity àtoms 

1/T 2 « (Au; 2 ) « C s , imp ((/io/47r) • 7l7s ^) 2 • ^ • ((S 2 ) - (S 2 f), (22) 

where 

(S 2 ) - (S 2 f « exp(-( 7s 7iI?/fcT)), (23) 

Cs.imp is concentration of magnètic àtoms in cm -3 , 7 s = 176radGHz/T, 7 i( 31 P) = 108radMHz/T, a is minimal 
distance of order of lattice constant, for Si a~ 3 w 5.0 • 10 22 cm -3 , 

(Alu 2 ) « 33.4 • 10~ 13 • (Cs,im P a 3 ) • exp(-( ls hB/kT)). (24) 

It follows that a direct suppression of decoherence due to interaction with magnètic impurity àtoms may be 
achieved by 

increase of relation (B/T) and 

decrease of impurity concentration Cs,i mp - 

To estimate the allowed concentration of paramagnetic centers in the silicon substrate Cs,i m p we will consider 
as previously the case slow frequency modulation and the decoherence time Td oc ~ ls. For B/T > 20 T/K we 
obtain 

1/T 2 « (A^ 2 ) « 0.74 • 10 3 • (C SMnp a 3 ) < 1 s~ 2 (25) 

or 

Csam P < 1-4 • 10" 3 • a- 3 w 0.7 • 10 20 cm" 3 . (26) 

Thus, the value of allowed concentration of magnètic impurity is practically unbounded at large enough value 
j s hB/kT. 

Another mechanism of individual qubit state decoherence is dipole-dipole interaction with nuclear spins 
I 7^ of impurity diamagnetic àtoms having concentration Ci^ mp . Isotope 29 Si with 7 i,i mp = -53 rad MHz/T 
is one of the such àtoms. 

In this case correlation function take form 

(A s ^(r)A Ws (0)) = (Alu 2 ) ■ exp(-í/T|| lmp ), r c = T||, irop , (27) 
where T\\ ^ mp ~ 10 4 s— impurity nuclear spin longitudinal relaxation time of isotope 29 P at low temperature |Q. 

47T 

For B = 2 T and for spin temperature Tj at which there are near-full polarization of impurity nuclear spins 

lUMpKB/ml > 1 (29) 

or for Tj < 0.8 mK, we will obtain that the allowed concentration of the isotope 29 Si in % is 

Ci,im P % < 4.5 • 10~ 2 % (30) 

This value may be increased provided a further decrease of spin temperature T\ . 
For comparison, natural abundance of isotope 29 Si in natural silicon is 4.7%. 

Hence the main reasons for the internal decoherence of qubit states are the modulation of resonance qubit 
frequency due to hyperfine interaction with fluctuating electron spin and due to interaction with randomly 
distributed impurity diamagnetic àtoms containing nuclear spins. 



1/T 2 » C I , imp ((/x /47r) • 7 i7Um P ft) 2 ■ ■ (1 - tznh 2 (\ lhimp \hB / m))) (28) 



5 The random phase errors in ensemble quantum register 



The preparation of both the initialized and other bàsic states of quantum register is followed by error generation 
through interaction of qubits with the environment. Any one-qubit error in either of the bàsic state is defined 
as a superposition bàsic states |0) and |1). 

As an example of one-qubit error action we consider the following unitary (it is not necessary unitary) 
transformation 



U K = 



1 



- e 2 + e 2 



1 + ie z 



1 - íe 7 . 



where e x , e y , e z are independent random error function of time. 

The action of this transformation on initialized state p(0) = |0) (0| = 
statc 



1 




p e (Q)=U e \0) <0| U- 1 = 1/2 
where the components of Bloch vector are 



P-.e 
1-P*. 



P B , e = l-2p e (0), P±, E = 2 



2 -4) 



exp(±i(^ £ 



(31) 



leads to the modified pure 



(32) 



(33) 



p e (0) = (e 2 x + e 2 y )/(l+el + e 2 y + e 2 z ), tan^ £ = (e x e y + e z )/(e x e z - e y ). (34) 

The fidelity of initial state |0) is determined by 

F(0)= Spp e (0)p(0) = l-p e (0), (35) 

where p e (0) is the random error probability. The diagonal elements of the perturbed density matrix p £ (0) 
depend only on the value of p e (0) and have no random phase factor. 

The random phase factors with arbitrary continuous random phase (p e have only non-diagonal elements of 
density matrix p £ (0) or the Bloch components P± e · 

All 2 L bàsic state vectors of quantum register |n) are generated by one-qubit unitary transformations iNOT 
which transform the initial basis states of certain qubits |0) in state |1). The containing error iNOT e operations 
leads also to appropriate random phase factors. 

The existence of the such uncontrollable phase factors can involve a problem, which mainly lies, as was 
noted by S.Kak Jn|, in the inability of the quantum error correction codes (QECC) to correct this analog-type 
quantum phase errors as "they can potentially correct only bit flips and phase flips and some combination 
thereof, which errors represent a small subset of all the error that can corrupt a quantum state" . 

We will advance here some arguments in support of the ensemble approach use in connection with the 
generation of uncontrollable phase factors. 

The producing of random phase factors (known as fasors) may be considered as a stationary random process 
realized by random functions s K (t), e y (t) and e z (t) possessed by a statistical ensemble. A time-average of density 
matrix is practically unfeasible as it requires prolonged state measurements of an individual quantum register. 
But if it is assumed that random process in ensemble quantum registers is ergodic random process, it may 
be suggested to substitute the time average by the ensemble average. Practically the ensemble averaging will 
occur by natural way, when unitary transformations and ensemble quantum register states measurements are 
produced at the same time for the whole ensemble by the same controlling pulse, much as in bulk-cnscmblc 
liquid quantum computer prototype. The ensemble-averaged reduced density matrix will describe now same 
mixed state whose non-diagonal matrix elements will decay due to decoherence process and have no analog-type 
random vàlues. Such states involve no special problems for QECC. 

The decay of non-diagonal matrix elements comes much rapidly the more dimension of non-diagonal matrix 
element blocks. Thus during the process of preparing of many-qubit quantum register basis states a ràpid 
decoherence can play the constructive role. The ràpid decoherence may be produced, among other processes, 
by means of a sharp decrease, for the short time, of external magnètic field. 

As the result mixing degree of reduced density matrix (p e ) may be treat as temperature analog in a liquid 
quantum computers. For suited conditions the value (p £ ) -C 1. Analogously to the low temperature many-qubit 



mixed state the considered state may be transformed to pseudpure state, but in so doing the entanglement 
property of large-scaled state at low temperature will not be violate In this case we will have the fully 

quantum computing instead of it imitation. 

I am grateful to K. A. Valiev for stimulating discussions and an inestimable support and to V. A. Kokin for 
the large technical hclp. 



Summary 

The main mechanisms of decoherence for low temperature are the adiabatic processes of random modulation 
of qubit resonance frequency without spin flips produced by hyperfine interaction of nuclear spin with clectron 
spin of the bàsic àtoms and dipole-dipole interaction with nuclear spins of impurity diamagnctic àtoms. It is 
estimated the allowed concentration of nuclear spin containing isotopes 29 P. It was shown that the random 
phase error of quantum register state in the ensemble approach are averaged and the pure bàsic state transforms 
to the mixed state, analogical to mixed state at nonzero temperatures. 
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